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Mueller Navelet jets were proposed 25 years ago as a decisive test of BFKL dynamics at hadron 
colliders. We here present the first next-to-leading BFKL study of the cross section and azimuthal 
decorrelation of these jets. This includes both next-to-leading corrections to the Green's function 
and next-to-leading corrections to the Mueller Navelet vertices. The obtained results for standard 
observables proposed for studies of Mueller Navelet jets show that both sources of corrections 
are of equal and big importance for final magnitude and behavior of observables, in particular for 
the LHC kinematics investigated here in detail. Our analysis reveals that the observables obtained 
within the complete next-to-leading order BFKL framework of the present work are quite similar 
to the same observables obtained within next-to-leading logarithm DGLAP type treatment. There 
is still a noticeable difference in both treatments for the ratio of the azimuthal angular moments 
(cos20)/(cos0). 
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1. Introduction 

The understanding of the high energy regime of QCD is one of the key questions of particle 
physics. In the semi-hard regime of a scattering process in which s S> —t, logarithms of the type 
[oi s ln(*/|? | )]" have to be resummed, giving the leading logarithmic (LL) Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) Pomeron [jl|] contribution to the gluon Green's function. The question of testing 
such effects experimentally then appeared, and various tests have been proposed in inclusive [Q], 
semi-inclusive [Q] and exclusive processes PJ]. The basic idea is to select specific observables 
which reduce the importance of usual collinear logarithmic effects a la DGLAP []|] with respect to 
the BFKL one: the involved transverse scales should thus be of similar order of magnitude. We here 
consider the Mueller Navelet (MN) jets ^ in hadron-hadron colliders, defined as being separated 
by a large relative rapidity, while having two similar transverse energies. In a DGLAP scenario, an 
almost back-to-back emission is expected, while the allowed BFKL emission of partons between 
these two jets leads in principle to a larger cross-section, with a reduced azimuthal correlation 
between them. We report on recent results where both the NLL Green's function [^] and the NLL 
result for the jet vertices $Q are taken into account 1 . These new results have been obtained based 



on a fast Fortran code, which allowed us to go beyond the studies of Ref. [ 10] where we developped 



an exploratory Mathematica code. Detailed results will be presented elsewhere |[T 
2. NLL calculation 

The two hadrons collide at a center of mass energy \fs producing two very forward jets, whose 
transverse momenta are labeled by Euclidean two dimensional vectors k/ i and k/ 2, and by their 
azimuthal angles tyj i and <pj2- The jet rapidities yj\ and yj 2 are related to the longitudinal mo- 

s/s 

da r 1 . r 1 N da a) 



mentum fractions of the jets via xj = ^e yj . For large xj\ and Xj t 2, collinear factorization leads to 



d|k 7a |d|k 7 , 2 |dv7i dy/,2 



V [ dXi / d*2/a(*l)/b(*2)-7ji 1 1 1. ^ , (2-1) 

tb^° J° d|k 7il |d|k 7 ,2|dy 7i idy /i2 



where f a b are the parton distribution functions (PDFs) of a parton a (b) in the according proton. 
The resummation of logarithmically enhanced contributions are included through kj -factorization: 

dOab 

dlkj,! ] d]kj- i2 ] dy^i dyj ;2 

where the BFKL Green's function G depends on s = x\X2S. The jet vertices V a ,b were calculated at 
NLL order in Ref. Combining the PDFs with the jet vertices one writes 

da 

d|ky ; i I d|ky j2 1 dy J; i dy J>2 

where <£>(kjj,xjj,ki) = J dXif(xi)V(ki,Xi). (2.3) 

In view of the azimuthal decorrelation we want to investigate, we define the coefficients 

^ m =j d0 7;1 d<^ 2 cos - (j) Ja - n)) 7d 2 ki d 2 k 2 <J>(k 7i i,^ jl , -ki)G(ki,k 2 , ^(kj^Xj^M), 

'These vertices have been recently recomputed in Ref. [B], in a full BFKL approach. 



ftyj,ityj? |d 2 k 1 d 2 k 2 K(-ki^i)G(k 1 ,k2,f)V b (k2,x 2 ), (2.2) 



= J d0 /;1 d0 7i2 J d 2 kid 2 k 2 0(k 7; i,^ jl ,-k 1 )G(ki,k 2 ,f)<I>(k 7 ,2,^,2,k2). 
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from which one can easily obtain the differential cross section and the azimuthal decorrelation as 

dcr ^ 

"o and (cos(m<p)) = (cos (m((j)j ;1 - (j) J>2 -n))} = (2.4) 



d | kj j | d | k 7j2 1 djy, 1 % 
The important step is then to use the LL-BFKL eigenfunctions 

^, v (k 1 ) = -L(k2)' v ^ e '^, (2.5) 

7l\/2 

although they strictly speaking do not diagonalize the NLL BFKL kernel. In the LL approximation, 

r ( s \ £0 ( m,v ) 

^ m = (4-35 m , ) J dvC m , v (|k/,ll»*Ai)C,v(|k/,2|,^)(-J , (2-6) 

where C my (\kj\,xj) =J d0 7 d 2 kdx/(x)V(k,x)£ m;V (k)cos(m0 7 ) , (2.7) 
and ®(n,v)=N c a s /KXo{\n\, 1 2 + iv), with 7) = 2¥(1) - V (7+ f) - ¥ (l - 7+ j) • The 



master formulae of the LL calculation (2.6, L2/7J) will also be used for the NLL calculation, the 



eigenvalue now turning to an operator containing a v derivative []12[ |13|], which acts on the impact 
factors and effectively leads to a contribution to the eigenvalue which depends on the impact factors. 

At NLL, the jet vertices are intimately dependent on the jet algorithm [§]. We here use the 
cone algorithm, with the cone parameter R = 0.5 2 . At NLL, one should also pay attention to the 

choice of scale sq. We find the choice of scale so = v/Sol so 2 with soi = T-kj, rather natural, 

v ' ' ' x j,i ' 

since it does not depend on the momenta ki^ to be integrated out. Besides, the dependence with 
respect to so of the whole amplitude can be studied, when taking into account the fact that both 
the NLL BFKL Green's function and the vertex functions are so dependent. In order to study the 
effect of possible collinear improvement fll5|], we have, additionally, implemented for n = the 
scheme 3 of the first paper of Ref. [|TJ]. This is only required by the Green's function since we 
could show by a numerical study that the jet vertices are free of 7 poles and thus do not call for any 



collinear improvement. In practice, the use of Eqs. (2.6, |2.7|) leads to the possibility to calculate 



for a limited number of m the coefficients C mv as universal grids in v, instead of using a two- 



dimensional grid in k space. We use MSTW 2008 PDFs [16] and a two-loop strong coupling with 
a scale fi R = |k /i2 |. 

3. Results 

Fig. H (left) displays the cross-section as a function of the relative jet rapidity Y, for the LHC 
center of mass energy y/s = 7 TeV, for which most of LHC data are taken at the moment, while 
Fig. [I] (right) shows the relative variation of the cross-section with respect to MSTW 2008 when 



changing the PDFs according to Ref. [17]. This explicitely shows the dramatic effect of the NLL 



vertex corrections, of the same order as the one for the Green's function [13, 18]. 



2 A detailed study pl[j, based on the work of Ref. [ |l4[ where the jet vertices were computed in an approximated 
small R treatment, shows that the difference between an exact treatment and this approximation is small. 
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Figure 1: Left: Differential cross section in dependence on Y for |ky i| = |k/2| = 35GeV, at yfs = 7 TeV. 
Blue: pure LL result; Brown: pure NLL result; Green: combination of LL vertices with the collinear 
improved NLL Green's function; Red: full NLL vertices with the collinear improved NLL Green's function. 
Right: Relative variation of the cross-section in the full NLL approach, for various choices of PDFs with 
respect to MSTW 2008 ones. 
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Figure 2: Left: Azimuthal correlation in dependence on Y for |k/,i | = \kj,2 \ = 35 Ge V , at y/s = 7 TeV. Blue: 
pure LL result; Brown: pure NLL result; Green: combination of LL vertices with the collinear improved 
NLL Green's function; Red: full NLL vertices with the collinear improved NLL Green's function. Right: 
Relative variation of the azimuthal correlation in the full NLL approach, for various choices of PDFs with 
respect to MSTW 2008 ones. 



Fig. ^| (left) displays the azimuthal correlation as a function of the relative jet rapidity Y, for 
the LHC center of mass energy yfs = 7 TeV, while Fig. ^| (right) shows the relative variation of 
the cross-section with respect to MSTW 2008 when changing the PDFs using the sets of Ref. [JO]. 
The decorrelation based on our full NLL analysis is very small, similar to the one based on NLO 
DGLAP. Fig. H (right) shows that (cos ) is much less sensitive to the PDFs than the cross section. 

Detailed studies [10] have shown that the main source of uncertainties is due to the renor- 
malization scale Hr and to the energy scale yfsQ. This is particularly important for the azimuthal 
correlation, which, when including a collinear improved Green's function, may exceed 1 for small 
Mr = Mf- A noticeable difference can be expected between BFKL and DGLAP type of treatment 
for the ratio (cos2<p) / (cos <p) . The BFKL results are shown in Fig. ||. We refer to Ref. [ |l0| ] for a 
comparison with the DGLAP approach at y/s =14 TeV. 

In order to get more insight into the azimuthal correlation between jets, we have studied the 
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Figure 3: Left: Ratio of azimuthal correlations (cos 20) / (cos0) in dependence on Y for |ky 1 1 = |k/ 2| = 
35 GeV , at y/s = 7 Te V. Blue: pure LL result; Brown: pure NLL result; Green: combination of LL vertices 
with the collinear improved NLL Green's function. Right: Relative variation of 
approach when using other PDF sets than MSTW 2008. 



(cos20) 



in the full NLL 



A0 distribution, a quantity which is accessible at experiments like ATLAS and CMS. Computing 
(cos(«0)) up to large values of n gives access to the angular distribution, since 



^ = ^{l+2lcos(«0)(cos(n0)> 



(3.1) 



The results, for various values of Y , are shown in Fig. ^| Our full NLL treatment, when com- 
pared with a mixed NLL+LL approach, predicts less decorrelation for the same Y, and a slower 
decorrelation with increasing Y. 
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Figure 4: Azimuthal correlation. Left: NLL Green's function combined with LL vertices. Right: NLL 
Green's function combined with NLL vertices. 
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